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Suspension of Solid Particles in a Viscous Liquid 


Identical spheres 
Concentration c 


Long wavelength assumption 
{AMP} = tkesks, Ps 


Effective medium properties 
Viscous liquid 
= 2 Z 
= —-l@ = 
(u 7) BA Be = dey +> Hep 
: 3 
Solid particles } Hef 
Pp Peff 


Effective properties depend on frequency, concentration, and scattering properties. 
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Modelling 


\ka| 


( 


MST 


|K.d| <1 . . 
% Mei’s self consistent model 


- Effective bulk modulus 


Hydrodynamic 
yorouy - Effective mass density 


Self Consistent % Ament’s hydrodynamic model 


[Mei, PRB, 2007] 
[Ament, JASA, 1953] 
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Objectives 
To investigate effective properties of a viscous suspension: 


inalvtical exoressions for scattering coefficients/for a spherical 


wave order In visco-acoustic medium. 


particle for arbitrary partia 


> to ext nd the self-consistent core-shell effective del|through 
inclusion of the wave mode conversion and calculate effective dynamic 
properties. 


nvestigate the dynamic effect of particle geometry jon the effective 
mass density by extending a hydrodynamic model. 
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. Motivation 

. Objectives 

. Single-Sphere Scattering 

. Effective Properties with Spheres _[Mei, PRB, 2007] 

. Effective Properties with Spheroids [ameri, sasa, 1953] 
. Conclusion & Summary 

. Perspectives 
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The Research Problem 
Why single-particle scattering? 


> Single particle theory 


> Multiple scattering theory 


> Effective medium theory 


Wave Mode Conversion 
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The Research Problem 


: Quadrupole Field 
> Lord Rayleigh (1871) » “ ~ Particle 
: : ? * 
= Rayleigh Scattering ka] <1 Incident Wave bs Forward Wave 
> Allegra and Hawley (1972) = \ ——- 
- Particle: solid; Host: viscous fluid (compression) : 
- Scattering coefficients order: monopole & dipole — monopole Field ape 


> Kuster and Toksoz (1974) 
- Particle: solid; Host: solid 


- Rayleigh scattering Ea ee 
- Scattering coefficients order: arbitrary 
> Pinfield and Forrester (2017) 


- System: same as [A&H,1972] (shear) 
- Scattering coefficients order: dipole 


> This work 
- Scattering coefficients order: arbitrary 
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General background 


Field equation for te )V(V-u)+ Wie Oru 
particle displacement, u Ht H —?P ar 


Helmholtz decomposition 


Scalar potential Vector potential 


Viscous fluid «1 = —iwn Time harmonic potential exp( — iwt) 


(V2+k2) = 0 mn eo a 
nee 1s \k.| < |k, 

“)>c) (V+kK)A=0 k=F? 
LL 


21 [Ke 
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Solution of the Helmholtz equations 
Spherical coordinate system (r, 0,0) r€[0,c0), 9€[0,x],o € [0,2 ) 


[V? +k? ]A =0 


Debye potentials 


(V? +k; )(vx) = 


Plane compressional wave © = ®, exp(ik,.r) = ®, exp(ik,rcos®) 


hl nee ee ea 


n=0 
® = yi (2n4+1)T"h, (k.r)P, (cos) a! = STi" (2n +1) Ae j, (kr), (cos) 
n=0 n=0 
= STi (Qn + i Bs h, (k.r)P, (cos) y= Sit (2n + 1) AS I (k’r) P, (cos@) 
n=0 n=0 
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Displacement components 


us = Va 
Boundary condition at r=a Boundary equations 


uf (r=a)=u;(r=al MX =D 
ug (x = a. | = Up. (x = a | M, My Miz Mig || tn D, 
A 7 My, My. My; May || Th" _ {Pe 
Ot (2 — a) ae (2 — a) Ms, M3. Msg Mg, |) At® Ds 
My, Myy My, My, |) Ass D, 


Q 

at 
——_ 
re 

| 

© 
VS 
| 

Q 

& | 
——_s 
re 

| 

© 
VS 
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Scattering coefficients 
Monopole n=0 Dipole n=1 


+h 
z. 


Quadrupole Field 
p - “| baa Particle 
/ 


Incident Wave 


* . 
enue? 


\ Senne? 


* es _— =~ 
Va mh ee 

/ 

! 

| 

\ 


Monopole Field 


Scattering profile of 
a single-sphere 
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Scattering coefficients 


oe 2i[2hy(y,)—ysho(y.)] 
i ee ree eae NIE ea 
135 3h, (y, ) + yt, (y, }| 


1 


Tss MS fee eee re Og ea D) 
> 93h, (y,) Fe (ye) Oy) 
2 SC n+l 
Ts = —-—______—_—______y? Te =O Ye 
: 3y, | 3h, (y, ) + y.b3 (y,)] Ye 
Ts _ O @) 
[3h, (y,) + yh) (y,)| 
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Numerical results 


. Yo Attenuation factor 
23. 10< 2.6x 10-19 
(Np.m-1.MHz~?) 
se Ai Frequency range 0.01 — 10 MHz 
a n+1 
jg (We) Ve fk,a] <1] |wetk,a) € {10-2, 10?) 


iG 


Sin (Ys) Re(k.a) = 100 > Re(k,a) = 0.04 
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Numerical results: Dipole coefficients 
-10° 
1047 
rs} 


10°! 


40°! 


0.03 
me = 0.02 
ow 
E 
(or) 
0.01 
0 
10° 10°! 10° 1 10° 10° 107 
Fte(k_a) 


10' 10° 
Kte(k_a) 


10° 
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Numerical results: Dipole coefficients 


RoC 


siaannnied Re(TS- me 


0.1 
0.05 
_ 0 
oz 
e 0.05 
= - 
67] 
-0.1 
-0.15 
HB 
10" 10°! 10° 10! 10° 40° 10°! 10° 10! 10° 
Fte(k_a) Rte(k_a) 
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Numerical results: Quadrupole coefficients 


CC 


107 10°! 10° 10! 10° 107 10°" 10° 10! 10° 
Rte(k_a) Re(k_a) 
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Numerical results: Quadrupole coefficients 


10° 
107 10°! 
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Numerical results: general results 


40° 


1 920 a 


404° L 


100° |) El Num. CS 
“&- Num. SC 
49°89 | | at Num. SS 
+ Analytical 


497100 
12345 6 7 8 9 101112131415 


Orders n 
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> Analytical approximations have been derived for scattering 
coefficients for general partial wave orders for a sphere 
In a viscous liquid. 


> Good agreement between the analytical and numerical solutions 
within the validity region. 
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. Motivation 

. Objectives 

. Single-Sphere Scattering 

. Effective Properties with Spheres 

. Effective Properties with Spheroids 
. Conclusion & Summary 

. Perspectives 
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Effective medium 


SS 


Effective medium 


%, Core-shell model: 
host medium 


> Concentration: 


particle 
= 3 
> Self-consistency: effective medium 
i — 0 
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effective medium 


Boundary condition at r = b Boundary equations 


My, My My Miy|] to R, 
ug (x =~ b) = Ug (x = b) My, My. My; Ma, |] ty | | Ro 
Ms, Msg. Ms3 Mg, || Ci Rs 
My My My My |] os Ry 
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Effective bulk modulus 


Monopole mode n = 0 Polar displacement and 
No shear Mave stress boundary equations vanish 
Viscosity in complex » and 


Solve radial displacement 
IX, and stress equations 
Ix.|<1 


|Self-consistency ti =0 | |Self-consistency ti =0 | = 0 


eS 


Te ox (k,b)" 
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Effective bulk modulus 


Kuster-Toksoz ca ead — 
3Bogp +4 2 3B + 4p Solid in solid 


l—c 
memes) a= 4 


Compressibility « = 1 / B Kore = (L—e)K + cK Volume-averaged 


[Kuster, Geophys, 2002] 
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Effective dynamic mass density 


> Rearrangement of boundary equations 
> Leading-order in X. and x, 


<1 , ; 
<1 Poff 


Dipole mode n=‘ 


2 
= 3( Tp — TTS )[x, bf(x,) - My ,)] 30" [x, 66,) — Zio) - 3x, 0 


X, 
+ iit [hy (x) x, mi&,)] +e |) —%, 1) 


2 
F = 6(TeTS — TPT" )[x, bi(s,) + 2b, (x,)] + OT" [x i{6s,) + 2i0,)] —3«,T" 6S 


>.< 
+ixtTS [hy (x,) — x, hi(x,)] +ix*[ I) —x, HO) 
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Effective dynamic mass density: low |k.a| form 


Te = xiey3 (6-1) a = ciey?y, (6-1) 


TS = icy, y? (61) They? (4 
9 : 9 


Poff ~ Agreement with the 
lim t8#£ =(1— 
ae 0 ( c) =< Ament’s static limit 


x, 0 


[Ament, JASA, 2005] 
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ka 


> 1 


Asymptotic form of _—— (1 4 i)Ba 
scattering coefficients : 
8 = Jpw / (2n,) > 1 


cc ip-l 3 
ale 0 
320 + 1 


Te ~ —ie s 20 + 1 ve ‘ 2 fOr ee 
/ / 
Ape ace 8) Tanenreromu 


~iy, A Peff __ 
—— fy (p + 2p’) — 2c(p' — p 


TS ~ -ie sin(y, ) 
[Ament, JASA, 2005] 


Page 28 /52 


Md Mahbub ALAM, PhD Defense, 12/09/2019 


<n 


| LE HAVRE 
Loughborough q 
University 


Effective density with spheres gy AO 7 


Laboratoire Ondes 
Milieux Complexes 


Numerical results 


kb] <1 Validity domain of the model 


kb 


= 


kb 


max 


250 


200 


Frequency/MHz 


0 0.2 0.4 0.6 0.8 4 
Particle size/micron 
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Numerical results: Variation with frequency at different particle sizes 


Resonance-like behaviour like MST 


e50nm 
4100 nm 


Re(P.¢/P) 


850nm 
4100 nm 
= 500 nm 
@1 micron 


@5 micron 
© 10 micron 
© 50 micron 


C= (0A 
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Numerical results: Variation with frequency at various concentrations 


< 
© 
@ 
a. 
— 
ic) 
co 


c = 0.1 to 0.6 
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Numerical results: variation with concentration at fixed Re(k,a) 


e**Re(k,a) =1 ee Re(k.a) =1 
— Re(k,a) =2 i w= Re(k,a) = 2 
=» Re(k,a) =5 =~ Re(k,a) =5 
=e Re(k,a) = 20 = Re(k,a) = 20 


Re(P.+/p) 


5 aia 
0.4 
concentration c concentration c 
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Numerical results: 
|Water_|Hexadecane | Sunfioweroil | Olive oil 


Sound velocity (m.s~‘*) 1497 1299 1470 1464 


Density (kg.m-?) 997 1000 920.6 915.8 


Shear viscosity (Pa.s) 0.000891 0.00663 0.054 0.092 


Attenuation coefficient 
0.023 0.145 1.15 1.625 
(Np.m—1.MHz-~?) 


Attenuation exponent p 2 2 1.77 2 


—— Water (viscous) 
| —-—Hexadecane 
aie Sunflower oil 

r = =Olive oil 


—Water(viscous)* . 
=== Hexadecane 


10° 
Frequency/MHz 


10° 10° 
Frequency/MHz 
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> Extension of a core-shell effective medium model by taking the 
viscosity of the host fluid into account through wave mode conversion 


> Calculation of the effective bulk modulus and effective density 
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Does particle shape affect the effective density? 


Different shapes = Different drags 


Prolate spheroid: deg < 4po) Oblate spheroid: eq > poi 


Incident 
wave 
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Ament’s model: Effective density 
> Conservation of momentum 
> Conservation of momentum density 


> Drag force = buoyant force 


Lane 
V p 


D = 6mna 1+ Ba — ial 1 + 2 
bp 1 
57 P Qn 
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2c(1—c)(p’—p) W i 2c (1—c)(p’ —p)Z 


Pete Ps — W? 47 w2 4.72 
R (Por } S( Pus ) 
[_Statictimie | [__tnviceld Himit_| 
w+ 0, S(Pee) > 0 w— 00, S(Pee) > 0 


rere (Por) — 0, Ww) — OO, W = 2(p’—p)(1—c) + 3p 


W — CO 


ears Ecce 
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Ament’s model extension 


Valid for particle of 
arbitrary shape 


D— D (spheroid ), Vo V (spheroid ), Pore — Pose ( spheroid ] 
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Drag on a spheroid 
Drag on a prolate spheroid [Lai & Mockros, JFM, 1972] 
1 1 
Sem | = 
c , 82mmay ec? [pw . |4_ F ja Bl a 


= sl oe +3 ae 3 Teg? 


E a 
1{ 1 l+e 1 
= 3{5 +1 og te|_2 ; ; 
— (max {a j61%eq = (min {apo %eq J) 
Bol 
Q, (x) = Etog| = +2] —1 


Limiting case: sphere 
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Effective dynamic density with prolate spheroids 
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Case 1:Ament a.) — a; cs > 0 


Peftp —? Pett Ps — 


Md Mahbub ALAM, PhD Defense, 12/09/2019 =————————_————————____ Page 42/52 


Effective density with spheroids “|= us _ we Ci 


Loughborough Laboratoire Ondes 
University Milieux Complexes 


Case 2: Static w — 0, 1/B—- co 


l 2 
S 1— 
(B}s00 |W, +22 / W, (1/8 }00 W? +2? 


P eftp a. aon (1 — c)p aie cp 
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——— 


Case 3: Inviscid WwW — ~, 8 — ow 


Paty _ Ep[p te(e’—p)|+(1—c)(0’ —p) 
p E.p + (1—c)(p'—p) 
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Effective density with oblate spheroids 
Drag on an oblate spheroid 


1 1 
om 32nna*, e* ip . |4 2 +t} 2 32nna, low 
D, = 8mma, | + ion 3 ea) 77)+ 7 3 — 
K K, 1) 


i 1 3 2 Yn 
eo eee mee 
E Jal 2 


a a, ; 
= [Ps = gq, (x) =1-xcot™ (x) = —Q, (ix) faa -[zee(a] 


2 


Same drag form ———+ Same effective density form 
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Numerical results: comparison of Ament model & GSCM 
1.12 ; 0.02 
1.11 
0.015 
oe —GSCM ~e 
1.09 | oo State = 0.01 
& ee Inviscid C= Q.1 5 
1.08 } 
0.005 
1.07 + 
06 be ee ae 0 
io" 3-10 10° 10' {07 10° 107 
Re(k_a) 
1.5 0.06 
0.05} 
—~ 0.04} 
Q 
~: 
c= 0.4 0.03 
= 
7-969 
0.01 
0 
jo" = 40 10° 40" 107 10° 
Re(k_a) 
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0.015, 7 7 
- - - Spheres ra me 
= = =Oblate spheroids | / ‘ 
= = =Prolate spheroids if ‘ 
/ 
i] \ 
—~ 0.01} iu ‘ 
— ip oe” 
‘o U < 
= He i 4 
I \ 
0.005 | ii in 
40 . 
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Bed \ at % 
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0 a i ii siiiil 1 iis 
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x 
_ Kf . 
Sz 0.03 Hl \ 
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> Calcuation of the effective mass density for a suspension of both 


prolate spheriods and oblate spheriods by extending an hydrodynamic 
model from spherical particles to spheriodal ones. 
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> Analytical expressions for scattering coefficients for arbitrary partial wave 
orders for a sphere in a viscous fluid. 


> Extension of the core-shell effective medium through inclusion of the fluid 
viscosity and determination of effective properties of using the model. 


> Investigation of the effective dynamic density with spheroidal particles by 
extending a hydrodynamic model. 
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> Calculating scattering coefficients for a spheroid taking 
wave mode conversion into consideration 


> Calculating scattering coefficients for a core-shell 
particle. 


> Calculating the effective shear/compressional 
wavenumber using the core-shell effective medium 


model 
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Why be interested in effective properties ? 


%, Determine the concentration of particles. 
For example, the encapsulated medics in an organ. 


% Predict average propagation in a complex medium 
and interaction with a defect as a tumor in an organ. 


+ Solve interface / Scattering problems 
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Effective medium 
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Effective medium 


Core-shell model 


host medium 


° _— particle 
‘ 


effective medium 


: a 
Concentration c = a 
b 


Self-consistency condition t?" = 0 
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Effective density with spheroids 
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